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bt Introduction 


In studies of the propagation of finite amplitude sound 


waves in liquids, recourse is generally made to the so-called non- 


eR aa ae eT 


‘linear parameter B/A, which can be Gefined by an expansion of 


ec nn 


: Ae 
eee changes in terms of density changes in the liquid for an | 
isentropic process. If we use a Taylor a _ for this relation 


2 
(P=Po) 
P= Pot GP o-p) + & on 2) le 
OyS Oy 


' 
; 
where P»P, » are the instantaneous and hydrostatic pressures, Bape 
Bic instantaneous and equilibrium densities, where the derivatives 


' 
lare taken at P = P, and constant entropy, and limit ourselves to 


terms through quadratic, we can rewrite the equation as 





? 
oO 








- B 
p= pe ra 5 * ate? ° (2) 


Since the first derivatives is the square of the infinitesimal 


sound velocity, oe » A= Poe » and we can write the ratio 
Be 2p) " (3) 
A 6@ ape 
p 
O O,S 


There are three principal method of determining B/A;3 
1) Finite Amplitude Determination, 
2) Static Determination, and 


3) Thermodynamic Determination. 


[ein te Amplitude Determination 





iIhis method makes use of the harmonic distortion ina 


a Dear 





tpropagatea wave. The amplitude of the second harmonic is compared 


°° tne awpiicude of the first harmonic. There are several 














Se] 
2 





different experimental approaches to the problem. ! 
iShklovskaya-Kordil used an acoustic filter to filter out the 3 
| 


second harmonic ana a quartz receiver to obtain data. Adler and 
i 
| 
| 


| 


mental techniques are discussed by L.K. Zarembo and V.A. Krasil 'nikowe 


Hi edemann> used the same arrangement except that the measurements 
were taken by determining the variation of the index of refraction 


of the liquid instead of using a receiver. Various other experi- 


mome problems associated with the above method are that 
‘the filter employed may interact with the system, and if it passes 
mone harmonic waves then it passes the fourth nemnonte also. 
Reflections may also be a problem especially if an acoustic 


‘receiver is used. Theoretically, the equations concerned with 


Bo ae 


‘finite amplitude are difficult to manipulate and, in some cases, 





‘cannot be solved without approximations on the non-linear effects. 





wtatic Determination 


A static determination of B/A requires precise measure- 


ents of density as a function of pressure at constant temperature. 











‘The data is analyzed for the slope and derivative of the slope of 
specific volume with respect to pressure. This information will 


k 


provide the non-linear parameter of the isothermal equation of state. 


parameter has been investigated by Coppens.* While the difference 
is quite small, accurate measurements of the thermodynamic 
‘parameters needed to calculate this difference are almost non- 


| 
The difference between this parameter and the isentropic non-linear 
| 
| 
‘existent. | ! 

4 


ee a ee 


Les tee ee Me Shed ha 
3 
+ 


The main difficulty with this method of determining B/A 


iis that very accurate data are required as the two derivatives are 

















A 
labout equal and are subtracted from each other. 
. 


Thermodynamic Method 





| This method involves measuring the velocity of sound 


rand its derivatives with respect to pressure and temperature. 
‘Because or the simplicity of this measurement compared to the other 





:determinations, the thermodynamic method was chosen for use in this | 
: 
' research. 

i 
| There are two devices which measure the velocity of sound 


| 


with the accuracy desired, the velocimeter and the interferomever. 





i ane velocimeter was chosen,as the interferometer, because of its 
movavle parts, would require elaborate seals to contain the 
pressures required. 


To establish the thermodynamic method, we rewrite 


f 


lequation (3) 6 


2 
B Coop) - Poe on a 2p 


@p* o 62 9 o 


(28 = 


0&0 ‘9 
the differential dS can be written: 


ome (oD) aT + (<3 _ a 


sing Maxwell's relations, we recast the above equation in the form 


|  ‘TaS = Cpa? = 1 Dp 


or an isentropic process 


' 














, | 
| 

| Now | 

i 

- it io 

| a Sq eo (SD) p oa CS) 9, pa »P ) 

| ! 

! 2 AV ac - | 

(F980? * GSP p gS psd? * (as) 7,p58 | 

| i 

Ac BVT /&e 

’ = = leans = 

| ép’s “opt, S Cp (Sp os 

| OV 

| where = Ae = 


Py 


1 


Substituting the above expression in equation (3), we obtain 





| oe Se +e bc 
POO an 
= O LC . 
Po t= Ts ., (SP P,s (3a) 
= (B/A)' + (B/A)" 
! C 
P 
where Cc = =, 
1 p pV 





Thus we have anrvequation where the components can be experimentally 


determined. 


) 
| 





The assumption that sound transmission is isentropic has 


been investigated by Pome enen He showed that if the absorption 


| 


rof sound is small over a wavelength (i.e., a/k << 1) and 


| wx(1 = 4) | 
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where ais the absorption coefficient, k the wave number, 6/272 the 


frequency, « the coefficient of thermal conductivity, cy (c.,) the 


iheat capacity for constant volume (pressure), Pp the density, and ! 


c the velocity of sound. In most liquids the absorption coefficient 





a is sufficiently small to satisfy Gol'dberg's requirements and in 


‘ 
the present case experimental observation indicates that ais 


sufficiently small here also. 
C33 
wa(L = =) _s 
) a << 1 
pec 
for the liquids investigated. Coppens3° also shows that higher 
Dr dier terms and the cross terms of the complete Taylor expansion are 


ineglipible. 
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DS Experimental Equipment and Procedure 


A 


Aw Moathod of measurins sound velocity. 





Tne method used for the measurement of sound velocity in 


|this research was the velocimeter method of Greenspan and Tschiegg® 
as adapted by Coppens.? The sample was contained in a cylindrical 
tube each end of which was terminated by a 5 Mc/s quaYtz. | 
efystal for the transmission and receptiontof sound pulses. A 
Jmowledge of the repetition rate of the sound pulses is sufficient 
Eo Getermine the time required for a sound pulse to travel the 


¢ 
q 


_ of the velocimeter. When the length of the velocimeter is 


| 
| 
| 


- known the velocity can be determined. 
} The principle of the timing method can be understood by 


means of Fig. l. 
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Pulse 
transmitting receiving transmitting receiving 
erystal (a) t=0 crystal erystal (b) ts? cryecal 
Second 

| Pulse 
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- : == 
ee 


| souna pulses are introduced at t = 0 and t = 2T and are detected by | 
ba second transducer at the opposite end of the velocimeter. The 
' Sine T can be adjusted so that the two pulses will arrive at the 
‘receiving transducer at the same time (Fig. 1). The repetition 
Site was aajusted so that the two pulses add for a maximum 


amplitude. 


Looxing at Fig. 1, a and c, we see that the repetition 


er cee ee 
ee Sane 


rate is the reciprocal of the time required for the pulse to travel 


‘ 
& 
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twice the length of the velocimeter, so that ¢ = Thi = Ce 


“LENE Ls TTS 


B. Description of equipment 

The equipment required to measure the velocity of sound 
las a function of temperature and pressure divides naturally into 
) 


| The pressure system provides and regulates pressures up 


four systems. 


to 2000 psi which can accurately be measured to within 2 psi. 
The temperature control and measurement system provides 
1a temperature range from Fo 6a to 523° for the sample when the 


pressure vessel is placed inside a high temperature oven. The 





easurement system consists of a platinum resistance thermometer 
connected to a Mualler bridge. A spotlight galvanometer is used to 


balance the bridge. 


Cae) re oe ay Pe a ee 


The electronic system produces the electrical pulses that 
are changed into the sound pulses by a transducer, accurately 
eeaures tne repetition frequency of the pulses, and displays the 
received pulses on an oscilloscope to allow a correct superposition 


tonoe made. 


| _ 
, a mn tent ttt tlM  gi  ——  —f ———a ——  — 











The pressure vessel and components consist of the 


wen 


wars 


velocimeter and associated equipment, and a container to allow 


pressure to be put on the sample. 


Dae cue vcveton: 
iit Shah Ne ee eh 


A pressure line from the vessel leads to a Heise-Bourdon 
wtube gauge and to a dry nitrogen tank which provides pressures up 
l to 2200 psi The pressure in the vessel is adjusted by a standard 
}nigh pressure regulator. The system also has a provision for 


4 e ° e 
' 2llowing alr evacuation of the pressure vessel by a vacuum pump. 


(2. Tempersture control and measurement: 
The platinum resistance thermometer is placed in the 


receptacle provided in the base of the pressure vessel. This 





it position inside the metal tends to average the temperature 


‘fluctuations of the oven and also permits a temperature reading 


purer e Of Fabien Eeest ~ 


near the sample. The galvanometer used with the Mueller bridge 
is a high sensitivity spotlight type. It has an internal resistance| 


of 78 ohms, sensitivity of 0.005 wamp/mm, and a period of 4.3 sec. 


3. Electronic system: 

A sine wave audio-oscillator provides the signal which 
leventually (when appropriately modified) excites the quartz trans- 
ducer. The frequency (repetition rate) of the oscillator is 
accurately counted by an electronic counter. The sine wave passes 


.tnrough a band pass filter into special electronic pulse-shaping 





units, athe band pass filter reduces electrical line interference 


saat afrvects the dist: .y wi wie oStli.oscopa. The power supply for 


the special electronic eae units was a constant voltage ) 




















aE omer supply. Direct current output was regulated to better than 
Nae Over a yam@e of 200 to 325 voltge ge eolated AC yoltage of 
| Gas volts is @lsot#vallable for tube heater circuits. Ine wine 
| Wave is first led to a Schmidt trigger where a train of reetangular 
pulses is produced at the audio oscillator frequency. Tne snape 
| of tne various signals are shown in Fig. 2. The signal then goes 
ito diccking oscillator #1 which produces spikes at the audio 
oscillator frequency. These spikes go to a blocking oscillator 
counter. Tne input signal is then used to excite an RC network 
jwhich does not allow the device to pass any other input Signals 
for a certain time interval. The time interval can be changed by 
external adjustments. This signal is then led to a monostable 


| 
multivibrator. The output here consists of a series of rectangular 








‘pulses. The frequency of the audio oscillator is a multiple of 
this pulse frequency. This signal and the signal from blocking 
oscillator #1 are now led to a gate circuit. The signal from the 
lmultivibrator acts as a switch to allow a certain number of spikes 
(usually 2) from blocking oscillator #1 to go to the final blocking 
oscillator. The interval between these two spikes is 2T and is 
determined by the frequency of the audio oscillator. The final 


blocking oscillator shapes the signal again, and the signal is sent 





* 
Oe ea a el, 


to the transmitting transducer in the velocimeter. The signal from | 





A. 


the receiving transducer is then amplified in the wide band chain | 





’ 


amplitier and sent to the oscilloscope for display. The trigger 


i for the oscilloscope is taken from blocking oscillator #2. 
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[te Prasesure wassel and componente: 


a 











“2né@ main component of the pressure vessel is the veloci- | 
Mtr. Pe guartz am'y line in the velocimeter is manufactured 
crow Amersill II quartz (an optical quality fused quartz)! see | 


fot. gse c+ 1S made in the form of a hollow cylinder with an | 








Wouter diameter of about 1 3/16 in. and an inner diameter of 5/8 in. 


Mmaere is a nole + in. in diameter perpendicular to the drilled 





| 
portion wnicn is used for filling and emptying the delay line. | 
aE encs are polished whiie the other surfaces have a ground Pinisn.| 
The maximum and minimum lengths ae diametrically opposed and 
Ting dure 1.998995 and 1.998930 in., respectively. The maximun | 
deviations of the ends from perfect planes are 30 and 50 x 107° ahs | 
‘Tne length of the delay line is greatest toward its center (around | 
lthe bore). ‘The linear thermal expansion coefficient a for quartz 
mo 7.2 x 107’/°C--10% while the isothermal bulk compressibility Ky 
ms) 18.6 x 1078 /psi+7%. These data were supplied by the manufacturer, 
The transducers are 1 3/16 in. X-cut quartz discs with a 


3/3 in. concentric disc of gold over chromium plated on one side. 


The metallic sample provides the electrode on the other side. All 


—————— 


surfaces are polished and the fundamental frequency is designed to 
be approximately 5 Mc/sec. 

The pressure vessel itself and the metallic components 
were machined from austinitic annealed stainless steel (type 321, 
134-7). doe G¥fe coes nov fFeachi Wits Che metals used. 

A oressure inlet is providec in the sice of the vessel 


Bae ~ oe 2 
et ooh a 4. i 
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Cap with Gasket Gasket 
Transducer £Or LOT 


and Gasket Cushioning Clamping 
Delay Line Transducer 
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SsoL.i@ materiai in the bottom of the pressure vessel. 
| The quartz delay iine, cushioned by gaskets, is placed 


Males 


fs a metaibic sleeve (see Pigs 3)— Caps containing electrical 

Eoniections nold the t@ansducers tienviy against the delay line 
Bwith the aid of a gasket to optain ciamped boundary conditions. 
th a hole in its center is positioned over the opening in 


he side of the celay line so tnat the metal can run into the 


ecavity upon melting. Since bismuth expands upon freezing, the 


) 


asove equipment must be able to rotate so tnat the metal can ve 
| 


poured rom vhe openings in the delay line into a receptacle under- 





imeatn the above parts. ‘This assembly is attacned to the pressure 
: aS - 

| cover by electrical leads from the cover to the caps. 
ian annealed copper gasket is used to seal the pressure vessel. 


m 


ihe cover is then bolted to the pressure vessel. 
\) 


| 





oem List of equipment: 
Pressure vessel and components 
Composition gaskets 
Mir: Armstrong Cork Company 
Asbestosenitrile rubber composition An-&890 Accopac 
Pressure gauze 
Heise-Bourdon Tube gauge temperature compensated; 
rated to 2000 psi 
High pressure regulator 
Mfo> Linde Division of Union Carbide Corporation 
| nese pressuce Titeizes an@ Hiving 


7° => * 4 “ a =< 


oa oN 
- , . ~ ma _-— -_ -~ a= Qara on + aa ~ ~ 
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Hign temperature oven 
ivr: American instrument Company 
a 


Oe on Ta as 
amiicod ec 


emwperature Bench Mode Oven: 


MasGmiun tenperaturerPradieuwas 2 C 


4 og 
remperature regalation: 0,3 °C 


“J 


7 

Ca®aice Nos 4+-3640 | 
MewanuM resistance thermmermeter | 
: 

: 





ifr: Rosemount Engineering Company, Model 162D 


“fr: Gray instrument Company, Model S-1003 


Type 74-15 delaying sweep with precision delay 
ee? Dlue on unit 


Type 76-02 dual trace 


High sensitivity svotiight galvanometer | 
Mfr: Rubicon Division of Minneapolis-Honeywell 
Regulator Company, Model 3432 
Wide band chain ampiifier 
Mfr: Spencer, Kennedy Laboratories, Incorporated, 
Model #202D, Bandwidth 100 ke to 200 Mc/sec. 
with a voltage gain of 20 a@b + 1 db. 
Oscilloscope 
Mfrs; Fairchild Dumont, Type 766-H/P2B 
Jee in Unit 
etpll tee: to Semv/on. | 
Mi..0 OSB Ll: ae 


%. ™ 


Mass liew.cvs Packard Company, Model: 200 CD | 








es al 


nD | 





leas 


BIiGCce orc COunver 


‘ir: Heplevt® Packaré’ Company, Model 523° DR | 


Constant voltase porer Buppi7 
“nm: associated Specialties Company, Model 3 ! 


Guartz celay line } 





' Miacr: Vaipey Corporation | 

. 

: - ° 

Ce. FC. mation of sample ane’ esmerinental apvaratus 
ane indium used was a sampie tnat hac been used 


Morevicusly in measurements of B/A, but since it had been used in 
4 
;tne presence of a nitrogen atmospnere in a quartz and stainless 





\steel container, it was considered reusable. The original metal 
ad been obtained from the Indium Corporation of America with a 
a purity of 99.999%. The bismuth was obtained from 
=... 


mmener Chemical Company (Cat. No. B-319) with an estimated purity 
of 99.98%. Analysis of the metal as certified by the supplier 





listed the following metallic impurities: 





pes onic. (As) 0.0000% 
Copper (Cu) 0.002 % 
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Before any cata for the alloys were taken, sound velocity 


Velocity measurements were wade at four concentrations 


thereware two 


cram, 


2a 
sectic points (at 34%, 66%, and 48% Bi 


i y “ 


immed). Dawa were récordedvat these pointssandsat ag@7Z Bi 


aA ~ 


Pencerntratvion. Gied insorderetoverwec wor 


. ~ ae - i. da 
AC Levee DOlme WWES ST 


fe de 
“ne aDp: ROaCheCiet Be 


alioy's chaxageteristice tosthat of dtugkaane 


meme OASUremMenTS were also talen for bismuth. 


orn 


welore any measurements were mace, the experimental 


= 


first thoroughly cleanec witn dilute nitric acid, 


iWaver, benzene, and acetone rinses. 


Z 


Tne a irst “ex ver imenta | 


4 


sanple was prepared by cleaning 


1 


tne metal as avove except that the acid rinse was with a much 


mMeaxer solution because the surface of the metal Tendec to oxidize. 





ine metals were then weighed in order to determine and obtain the 
{ 


Wesired percentage alloy. They were then put into a heaker and 


T 
Ineated uncer paraffin oil to about 300°C, 
; 


Riz 
fenis point for about five 


THe ni xtureuwas held ay 





minutes and was stirred vigorously. The 
alloy was then allowed to cool while the stirring continued. The 
Wireezing point of the alioy was determined roughly and checked with 
line phase diagram. The alloy was then cleaned with the above 
aie tioned rinses and placed in the cup above the delay line. 

pince 2naium iS a fairiy higm priced metal, a technique 
Mapes ucCGssivée dicution was used. she samples were tested in order 


Sample had been 


cleaned and weighed for the 
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were obtained for indium which agreed within 0.6 m/sec. 
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Percer.t bismuth was 67.00%. 


a a A ane 
On Git? conc en er cl 


C 
| 

wL7@ Bi.), 
iC(hereafter designated as 48% Bi.), 


as 668 Bi.). 


33.84% Bi (hereafter 





inte OLGEer 


‘Sample, a wetting agent was used. 


aa 


Breuic (cs 350, 12,500) serve 


tela lS 
Vides 
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! 
The 
) 


transducers were thor 


ecetone rinses. 





wae coating had to be very thin in 


4 i + 
iwelcclty Cata. 
inte 





ference patterns (Newton's 


were tnen assembled and the whole assembly was placed in the oven. 


— 
fi 


ithe over 


vessel) 


| 


syeossel and 


cvne desired temperature. 
hewn OC AGms eel. 
ep. wa MA TPO EN 


Se Ee rs oe 
Oe a" BO 


—— aD wa ~ - 


Mer avith th 


migue describec apove 


ewe Ores 
sion values is considered to be about 1%. 
moutedé concentrations were 17.00% Bi ( 


Gesignated as 34% 


to gece the acoustical signal inuo and ou. cofec 
g 


oughly cleaned with acic, benzene, and 


They were then coated with the wetting agent. 


in all cases tne coating was tnin enough so that 


rings) were seen. 


+, pressure, anc thermometric connections were made. 


was brought slowly (to prevent warping of the pressure 


tase Wes come 
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After the sample had melted (as indicated by the 


a ‘ - 4. te = 


@ypearance of a transmitted signal on the oscilloscope), several 


sa S -~ oT} =~ ae we ~ A -* ~ ” ~ “ 7 st. Cs. 
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@nc 3 nours™ sotween temperature readings. t was found tiat efeee 


ition rate hed reached 90% of its finer value, | 


rene 
Batter 2 nours 99% of its final value; so that in a three nour 


Heal ©) ght nour perioag was ailowed to elapse for a temperature 
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Seei0d, tne equilibrium error was well Within r.m.s. reading error, | 


ichange. Data taxen for a combination temyerature and pressure i 


| 
ehange incicatec that in a four nour period the repetition rate 


mas 99% of the final rate. 





me readings were taken by adjusting the repetition 


mie. 1) was a maximum. This procedure insured that the time 
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irate so that the superposition of the two pulses at t = 3T (see 
Oecween che pulses was the time requirec for one pulse to travel | 








tO One end and return. The wave packet at t = T was also displayed 
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On the osciiloscope to aid in cetermining correct superposition. 


the frequency of the audio oscillator was then determined from the 
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Memes six days above melting temperature, a signal of correct form 


appeared.. Towever, the empiituce kent decreasing each day so that 
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‘III. Corrections and Calibratio 1 


The equipment's design provides a high degree of accuracy | 


0} 


DT Pa 


in measuring relative sound velocity in the delay line. However, 


—— 


s) 


he 


errors in the absolute determination of the velocity of sound may : 
rise from two different sources. The instruments such as the } 
chermometer, pressure gauge, and electronic counter should be | 


calibrated with the absolute values of their respective quantities 


om ee ee ee a 


SO that absolute values of the sound velocity can be measured as a 
| function of absolute quantities. The fact that measurements are 

Weng place ina Gelay line instead of free space also influences 
;these measurements. While some of these conditions are systematic 


respect to pressure and temperature, an attempt is made to correct 


awe. 


tor all possible errors. 





| A. Delay line correction 
| Since the length of the delay line changes with 


temperature and pressure, an error is introduced thereby in the 


= will not greatly affect the accuracy of the derivatives with 





determination of the velocity. It can be shown that the length 
9g 
change is governed by the equation 
L= Lt [ita(t-20°C) ]f. - + Kap] (4) 
3 | 


| 
| 
The maximum change in length of the delay line amounted to 0.0017 cm | 
which caused a change of 0.5 m/sec in velocity. Because of the | 
uncertainty in the v..ues of G@ and Km , the maximum uncertainty due 
| : 


‘to this correction was about 0.003 m/sec. 
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(a) Waveform of a pulse at t=T 
(b) Wavefore of the above pulse at t=3T 
Ihe experimental wave packet is shown by the continuous 


lines (/)\). 


The calculated wave packet is shown by the heavy dots | 
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Sirst half of the third cycles The only error intredue@d by this 


sare 


superimposed as shown. The only pulses that could be added are the | 
if 

t 

: 

| superposition is due to waveguide effects. However, since these 


pulses are smali, the changes in the sum of the amplitude is also 





| 

| mai. Better accuracy can be obtained in finding a maximum - 

| amplitude by superimposing pulses of larger magnitude. Therefore, 
| the pulses with the arrows drawn ahove were superimposed. A 


-correction to the velocity of sound must now be made since the 


-ccond pulse enters the delay line a certain time later in order 





ito superimpose at the new position. Since the delay is one wave- 


‘length over a distance of 24, the correction to the sound velocity 





its 

Ac = ce (5) 
Ac varied from 5.225 to 9.664 m/sec, with an accountable uncertainty 
iof 0.003 m/sec. One sound velocity measurement was taken using the 
superposition and compared to one using the above super- 
iposition. Differences in velocityamounted to 9.454 m/sec, all of 
which could be accounted for by equation (5 ) except for 0.210 m/secl 
Experimental uncertainties amounted to 0.177 m/sec, so the unaccoun 


lable uncertainty this "incorrect" superposition introduces amounts 


to about 0.04 m/sec. -This uncertainty is probahly due to waveguide 





‘effects. 


TE. 


C. Waveguide effect correction 
Waveguide effects account for the next largest 
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i each cycle does travel at the free field velocity, but since 


frequency. This produces a dispersive effect. The leading edge of 








measurements were taken at maximum amplitude, this effect will 


cause a systematic error. The amount of this error can be 


7 
| 
| 
| 
| 


e@culated for a particular arrangement./¢ The error in the 





“measured velocity of sound due to waveguide effects varied between 


10.160 and 0.380 m/sec. The uncertainty introduced by this correc- 





cion is negligible. 


D. Viscosity correction 


a 
o 
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i 

| 

: Viscous forces in a medium tend to slow the passage 
| sound waves. Kinsler and Frey state the expression for the 


velocity change Ac asl3 


| 
2 
| iee= es] (6) 
bap va? 
where u/p is kinematic viscosity, v is frequency of sound, and a is 


ithe radius of waveguide. This correction ranges from 0.100 to 





0.049 m/sec. Values of kinematic viscosity were found only for 
indiun? In BiL?, and bismuth.’ Values required for other concentra- 
tions were calculated from above values. Assigning an error of 50% 
sO the value of kinematic viscosity used gives an estimated 
iertainty of 0.004% m/sec here. Another possible source of error 
in the interpolated data may be the variation of viscosity with 
pressure,!* but since no data could be found on its magnitude in 
hiquia metals (in organic liquids, the viscosity roughly doubles 
ror a pressure change from 1 to 2000 Kg/em~) no estimate was made. 


Since the viscosity correction is small and a 50% error was assumed 


rE 
jin the value of kinematic viscosity, any variation due to pressure 
Dnata tment orn ees 
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tanould lie within these limits. 
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bh. Jemperature errors 








! The maximum temperature error is dependent on the 


jaccuracy or the calibration of the thermometer because the fluctua- 
; 
tions in the reading of the bridge were within the accuracy of the 


calioration of the platinum thermometer. 


| KX, Calibration 
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| In order to determine the absolute velocity of sound 


t 
chermometer, pressure gauge, and electronic counter was checked 


a re es es 


| 
as a function of temperature and pressure, the accuracy of the : 
‘against absolute values. 

The platinum resistance thermometer was calibrated by 
Rosemount Engineering Company and found to be accurate to within 
+#0.015°C from ~183°C to +350°%, 

The Mueller Bridge has a guaranteed accuracy of +0.005% 
+0.00005 ohms at a temperature of 23 eC. 

The Heise-Bourdon Gauge was calibrated with a dead weight 


piston gauge and found to be accurate within +2 psi. 


The electronic counter was calibrated with WWV and 


adjusted for correct readings. It has a stability of two counts 


; 


in one million per week and is accurate to +l count. 
| 
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3: Data 
| 


Ae Sapeoc colevlation for velacity "as so ° 


nm 





dA calculation is shown imedetail below Hor Gme cL Geae 


| points of the 34% Bi sample. The experimentally determined data 


Repetition rate 20.76668 +0.00077 kc. 


| 
| 
| 





2108.892 m/sec. 


Temperature 77 ot8 +0.02°C 

Pressure O + 2 psi | 

Actually the temperature is not determined directly, but is 

Fobtained from a table where resistance values are listed opposite 

cne temperatures. 

(Corrections for changes in the length of delay line: 

2t= 2k firo(t-20°C) If - > Kop J (+) 
= (10.154752 cm) (1.000041) (1) 

' = 10.155172 cn. | 

) 

) Hence 

) c = 2fy = (10.155172 cm) (20.76668) ke. 

| | 
4 
| 


Correction for superposition: 


ne = he oe - (2108892) 
2 2th (10.155172) (5x10°) 
= 8.759 m/sec. 
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Correction for waveguide effects: 


A 
B . (0.62) (x) 
“i am r 


——— oe 


x - Distance pulse travels 
a - Radius of delay line 


A - wavelength in medium 


(B/a) ae = 60.92) xe = 1.35107 )e = 0.922171 


(On6':)( 5°10") 


re ee a a 
eine ~—w 


1(8/w) - p/w)” 


| 9 = cott = 1.10003 radians 
1 + (B/u) 
2 
| 1 H(B/u)” + © 
: wty = a a i. 195793 
| 1 + y8(B/u) 
wtp = 1.428036 
4 
iG 
Cc = 758 =" 0.3 3am/cece 


mPorrection for viscosity: 


3 ss s 
ene —1.0-30°" _ = 0,042 m/sec. (6 
- elt 23 Tn atio (0.6h) m/ BOGS: AGe 


mp va 


(B/w) = 0.307390 9 = 1.41709 radians 
x=5 


Velocity of sound corrected for changes in length of delay line 


; 2108.892 m/sec 


‘Correction for superposition 8.759 m/sec 
i 
he 


Vf 
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mOormectio:n Cor waveralicde effect 0.331 m/sec 


0.042 m/sec 
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Brror due to uncertainties in temperature 


| 

| 

\ and pressure corrections of delay line: 0.003 m/sec. | 

| Error due to uncertainties in superposition | 
correction: 0.003 m/sec. 





| 
Error due to uncertaintLés in viscosity 


| 
correction: 0.004 m/sec. 


‘Srror due to root mean square Geviations in 





| obtaining repetition rates: 0.078 m/sec. 
feereeted sound velocity 2117.362 +0.088 m/sec. 


: 
, The experimentally determined sound velocities of the 
j274 Bi concentrations are given in Table I. The results of the 


alculations for each set of concentrations with the same pressures 
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were fitted by a least squares method to polynomial equations by 


. Soke 


the Brown University IBM 7070 computer. The results of the 





= for the sound velocity of the 17% Bi sample are also 


given in Table I below the experimental values of c. Polynomial 


ee Cel 


curves of degree one through three were calculated. The linear 
east square fit resulted in a smaller root mean Square deviation 


rom the experimental points of the 17% Bi alloy than the second 





and third degree polynomial equations while the second degree 

a least square fit resulted in smaller r.m.s. deviations 
from the experimental points of the 34% Bi, 48% Bi, and 66% Bi alloy! 
than the first and third degree one However, there was 


Still quite a large curve fitting error, even larger than experi- 


penta error in some cases. 





| 
4 Table II, wnich tabulates sound velocity as a function | 
bot temperature and pressure was computed from the polynomial 





Ose eee 
o-Lhel o um 


1 
EE SP er ——hr =r rr en | ae wea 








equations in the case of the alloys. 
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given are the experimental results. 
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For bismuth, the points 
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BD. Velocity of Sound Tables 


TABLE II 


Jalues of the velocity of sound (meters/sec.) as determined from 














22202 4:19 
Cen e525 
2212.630 
2207.464 
2201.482 


2117 485 
PlieeooG 
21000261. 
2104.573 
2100.727 


Pressure (Psi) 
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Pressure (Psi) 


a2 


Pemp 0 500 1000 1500 2000 

| 
| Bismuth | 
}31.8 16392433 1640.174%  1640.943  1641.866 1642.542 
328 1637.768  1638.552 1639.206 1639.927 1640,612 
342 1637.040 1637.647 1638.191 1638.991 
ape? 6330033 
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| Since the errors for all but the deviations when obtain-~ 
ling repetition rates are very small, they can be considered the 


‘same ror all points. The r.em.es. error varied from about 





0.030 m/sec to a maximum of 0.075 m/sec. The curve fitting error 
lalso varied widely from a negligible error up to 0.137 m/sec with 


i 
lan average of 0.52 m/sec. 





_— earlier discussion: 
Error due to uncertainties in 
temperature and pressure corrections: 0.003 m/sec. 
Error due to uncertainties in super- 
position corrections: 
Error due to uncertainties in viscosity 
correction: 0.004 m/sec. 
Error due to root mean square deviations 
in Be aaa repetition rates: 


Average 0.075 m/sec. Maximum 0.116 m/sec. 


0.003 m/sec. 
| 


‘Total error of experimental points 






Onl 26am cece 


Maximiim curve fitting error 0.137 m/sec. 
Averace curve fitting error “ m/sec. 
Serre enn ere mene ———_—. ao 
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Average 0.085 m/sec. iuaximum 
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C. Computation_of Be 





After the equations of sound velocity at constant 
® 
pressure were determined as a function of temperature, equations 
| for sound velocity at constant temperature were determined aS a 


§ 
De 


function of pressure. These curves were plotted and found to have 





a linear dependence. Graph V shows the curve with the maximum 


et eee 


Slope and Graph VI the curve with minimum slope. 
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Beside the value of c and its derivatives, the computation 


a er 
Te 


of B/A required three experimentally determined values. These are 


censity, thermal expansion coefficient, and heat capacity of the 





samples. 

| Values for the density of liquid In-Bi mixtures could not 
Ibe found. <A density curve vs. %Bi(by weight) at room temperature 
eur for the solia)t® showed an approximately linear relation. 

qua tions to find density as a function of temperature? were 
extrapolated and density values computed using a linear relation. 
However , Since bismuth expands upon freezing and the freezing point 


of ali the alloys tests was below that of bismuth, an attempt was | 





made to determine the density and change of density as a function 1 


u 
Le ° 


of temperature for the 66%Bi sample. The data obtained agreed 













toughly with the calculated values. 
The thermal expansion coefficient was not found for the 
lloys, but since indium anc bismuth each have the same value for20 
| he coefficient and the term containing this quantity is by far the 
on aller of the two making up B/A (Eq.3a), negligible error should be 
Hntroduced here. 


lieat capacity (C,, ) Mages gor the @lloysialso were 


unavailable. A iinear wats verdpeion between the values 
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found for bismuth and indium was therefore usedJ/ Any error in this | 





| value should have a very Small efféct a6 At is also in thé Gmailér 


erm of B/A,. 


sne results of the caiculations for density and heat 


‘ CR ET, PE PS 


Loca pacity are listed in Table IV. Values of B/A were computed from | 
\ 
os Values iistead in Table IV according to equation (3), and are | 


tabulated ioe the Cable, 
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| The error in determining B/A for the alloys is mainly the 
ninth of the values of (B/A)" since (B/A)" has @n order of ) 
\ magnitude equal to the uncertainty of (B/A)'. Due to the estimated | 
eeror in the density valmes (+0.5 gm), an uncertainty of 0.15 is 
scalcuiatea for B/A. The uncertainties in velocity are very small 


he the maxémaum uncertainty introduced to B/A is 0.00003. The 


tne siope of c vs. pressure affects the accuracy of 


f 
© 
© 
3 
cr 
a) 
pte 
ae 
& 
b Jj 
=s 


iB/A greatly since the slope of this line is very small. Errors in 
the data points on the order of +0.1 m/sec. will cause an 
‘uncertainty of 0.3 in B/A. The dependence of c vs. p appears to be 


inear well within experimental error, so this assumption was used | 


bs 


n the computing uncertainties above. The use of the average value 


of r.em.s. errors and curve fitting errors resulted in an average 





E 
Eee ertainty to B/A of 0.218. The use of the maximum values 


fesulted in an uncertainty of 0.58. 


eR er vee 


The calculation for the uncertainty for bismuth follows 
whe same procedure as above, with a much smaller uncertainty due 
re Gensity and a much larger uncertainty cue to the slope of ¢c vs. 


S as Gata at only three temperature points were used. 


TE 


| Errors for alloys 

\Zeror in (B/A) due to uncertainties of -asity CaS 

t 

f 

\Srror tn (B/A) due to uncertainties of c 0.0008 
} 


[aor in «P/A) cue uncertainties of (Se/ép).(av) D206 (mse 


0.433 


See uncer<ainty in (3/4) (av) 0.37 (max.) 
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L-rons for Bissmiatn 


4 
WSrror in (3/4) due to uncertainties of density 0.09 
| Error i (o/e) Gue to uncertaintiecmoure 0.000% 
tf : . 
"Seror in (B/A) due to uncertainties of (Oc/dp), 0.96 
f 
§rotal uncertainty in (B/A) for bismuth 2 
, | 
= | 
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| 
ie Conclusiong | 
; As far as could be determined there had been no previous | 
{ 


work done in determining B/A for indium-bismuth alloys. B/A 


== we ee ee we - _ 


measurements for pure indium and velocity measurements for bismutn 
as a function of temperature have been reported. The available 


Se ts with the appropriate results of this work are given below. 











le 
\Temp (°C) 160 161 
(e/ ep) ,x107°(em/ayne-sec) 1.47 149 
(se/et), (cm/sec- °C) ~28.5+0.6 -27.4 10.2 
le (m/sec) 220, 22-066 2220 Bi 
CB/A)? 4.78 502 
ee" Or: -0.26 
Mpa) Le 545% 4.8 

Bismuth 
Temn (°C) Velocity (m/sec) (Qc/at) (m/sec°C) (Source) 
. ola Oso. -~O.5 (20) 
| 285 1663+5 (22) 
| 289 1666 0 (21) 
B10 POS ee -0.16 (This work) 


ane results of this work appear to approach the values 


sr indium as the concentration of indium is increased (see also 





Meecions 1V-3 and D). The velocity of bismuth found in this 


| v 
a ‘* P ” t.* 5 * . 
me use. Lo extrapolate to the temperature usec. Otherwise the 


es fe jier »-iueen the t.- “ances of reported values. | 


h 
nceesearch agrees with Kleppa's value if (8c/d%t)_ found in this work | 
| 




















i L% is possible tnat the sound velocity measurements 


mS 


eported Kere were made under non-equilibrium conditions or condi-=- | 

tions wnich were now considered when malimag applicable correcvions. 

ror about six Gays after the sample had melted, the signal on the . 

hoscilloscope indicated poor acoustical contact of the sample with 
| 

the transducers. inen the signal changed to one indicating good | 

acoustical contact. However, the voltage amplitude of the signal | 

Searted to decrease slowly, so that after five more days an 

accurate superposition could no longer be made. Upon disassembly 

of the pressure vessei and components, the transducers were found 


+ 


vO be cracxed. Bismuth hac leaked out of the delay line and was 


| 
| 


found in the sleeve between the transducer and cap at one end of 


WO 








ithe velocimeter. The reason for loss of amplitude of the signal | 





was probably due to this leakage of bismuth. The bismuth collect 





ing behind tne transducers may have affected some of the earlier 


nw ee 


assumptions used in computing corrections such as the assumptions 





lof air-backxed transducers, a filled delay line, and a fundamental 
frecuency of 5 Mc/s. The bismuth probably leaked 

through a crack in one Of the transducers caused by thermal 

Meresses set up by the high temperatures. Due to these uncertainties 
n error in velocity measurements of 2 m/s is assigned. A larger 
error is not considered reasonable because the waveform observed 


'did agree very closely with the calculated waveform (see Section 





jiIi-3B), bus since the set of velocity measurements was not completely 


“9 


vhe avove esvimave or error was made. | 
nw aorlavtsevy derived an expression for the velocity of 


| 
moma, in Wexcures as u. .wmcrion of concentration of the components.* 
} i a el, 

















—— ee 
if 

{| 

| - 
SAssuming an ideal mixture, we have 


(7) | 


1“ is tne molar fracction Of age 


Bwnere U is vne internal energy and 
Beaifst component. A#1so 


| V = hiv, + (1-1)v, 


_wnere V is the volume. Differentiating equation(7) twice with 


}respect to volume at constant entropy and noting that 


| any oi 
eve § ve 


iwnere M is tne molecular weignt and c is sound velocity, 


r 


) My 9 io 
c* a see, + (eee ce 
M se M 


where M is the average molecular weight. A plot of the above 


ivalues of sound velocity is showm on Graph 7 along with the experi- 


‘ 


menval results of Coppens and tnis * Seal The maximum difference 
between values of c is 3%. 
| Tne reason the c vs. temperature dependence for three of 


‘ 


Ene alloys has a varying slope is not known but may be cue to the 





.aifference in slope for this dependence in the pure metals or may 


! 


iincicate a change in tne siope of temperature dependence of bismuth 


Y 
8 
i : 
Mbeyond the range investigatec. 


17 


ine B/A data for %he indiumebismuth system are shown on 


Sctepn &, The values increase regularly with concentration with the 


emma arte 
' 





jeoxception of the vm BiaAlloy, ‘The compound In.Bi is formed, and 
g a 


eee: Mek Come by Depettioin this apeg of concentration. t 
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-No other unusual characteristics have been reported. 


The non-linear parameter of the equation of state | 


‘even a rough guide been established for predicting values. Thus 


;far, experimental values have been obtained for some metallic 


iF = p(p,s) has not been derived from first principles, nor has | 





selements, alloys, water, and various organic fluids. Reported | 


k 


: 


' 
; 


femies of B/A rasge from 2.7 to about 1005. 


’ 


| 
There are a few physical constants that appear to be : 
| 





related to B/A. The metals provide a range of B/A values from 2.7 


jto 7.8. Physical constants that follow this trend are l/c, 


> 





IRL 


il/latent heat of vaporization, and density (see Graph 9). However, 


4 


h 





constant to follow the trend was 1/c. The fact that l/c appears 
ito be proportional to B/A is somewhat surprising as in equation 3a, 
it is directly proportional to B/A. Acoustic impedance (pc), ao 
Bnd B do not follow this trend. Table V lists some of these 
parameters. 

The effect of temperature on B/A may be important in 
determining what physical mechanisms determine the degree of non- 


linearity. Water changes by a factor of 1.9 over the temperature 


range GOO by while ethanol's maximum variation in B/A is 0.6. 





The fact that there are less parameters that follow the 
‘ 


when water and organic liquids were considered, the only physical | 
f 

| 
trend of B/A for organic liquids is not considered unusual, because 
| 





the organic molecules are much more complex so that combinations of 


m 








| 





‘Material 





$ 


Methanol 
peehy- 
Acetate 


‘Acetone 
5 
g 
3 


Diethy= 
‘lamine 


n-Propanol 


'Ethanol 


; 





Ref. 


Ref. 
Ref. 
Ref. 


Ref. 


Temp 


(°c) 


161, 
24.0 
160 
10 
318 
30 
30 
20 


30 


O 


7 
ng 
25 
26 
27 


B/A 


ae 
gy? 
4 50 
46° 
Pal 

7.8" 
9.04% 
eee 
9. 5° 
9.74 


10,3¢ 


10,52 
10.54 


TABLE V 
Physical Constants of Some Substances 





Velocity Latent 


(m/sec) eae of 
tion 
(cal/gm? 
ele 10058 
24-70" 573° 
2313» 4.682 
14.9° 5908 
1639 oon 
144.0% Ton 
12754 94 @80° 8 
1188¢  125@56°8 
1189'  263@65° 8 
11287 =. 124@58° 
1095¢ 91@58°8 
1295  164@97° & 
1232! 204078 & 
f Rete 26 
Ref. 29 
h, Ref. 20 
7 (het. 13 


Density 
( gm/em>) 


OBE 
6.95% 
7.03% 
1.008 


Og 
13552 


0.888 
0.798 
0.798 
O. 928 


oO. 716 


0.808 
0.798 


esl 


Isother4 
mal 
Compres- 
sibility, 


(em? feyre) 





o1.08 | 
3.1 on 
2nee | 
| 

al 


10225 
12.88 


1O.7= 


10.58 
9.98 
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